setup;
global h

x = 2 : h : 8;
y = 2 : h : 8;

m = length(x);
n = length(y);

u = ones(m,n);
f = zeros(m,n);

g = zeros(m,n);
g(1,:) = u(1,:);
g(:,1) = u(:,1);
g(m,:) = u(m,:);
g(:,n) = u(:,n);

tic
v = solverPoi(crop(f), g);
toc
hold on;
figure; surf(crop(u) - v);
figure; surf(v);

abs(crop(u) - v)
